HIGHER-ORDER GLOBAL REGULARITY OF AN INVISCID 
VOIGT-REGULARIZATION OF THE THREE-DIMENSIONAL 
INVISCID RESISTIVE MAGNETOHYDRODYNAMIC 

EQUATIONS 



Abstract. We prove existence, uniqueness, and higher-order global regularity 
of strong solutions to a particular Voigt-regularization of the three-dimensional 
inviscid resistive Magnetohydrodynamic (MHD) equations. Specifically, the 
coupling of a resistive magnetic field to the Euler-Voigt model is introduced 
to form an inviscid rogularization of the inviscid resistive MHD system. The 
results hold in both the whole space K.^ and in the context of periodic bound- 
ary conditions. Weak solutions for this regularized model are also considered, 
and proven to exist globally in time, but the question of uniqueness for weak 
solutions is still open. Since the main purpose of this line of research is to 
introduce a reliable and stable inviscid numerical regularization of the under- 
lying model we, in particular, show that the solutions of the Voigt regularized 
system converge, as the regularization parameter « — >■ 0, to strong solutions 
of the original inviscid resistive MHD, on the corresponding time interval of 
existence of the latter. Moreover, we also establish a new criterion for blow-up 
of solutions to the original MHD system inspired by this Voigt regularization. 
This type of regularization, and the corresponding results, arc valid for, and 
can also be applied to, a wide class of hydrodynamic models. 



with appropriate boundary and initial conditions, discussed below. Here, v > is 
the fluid viscosity, /i > is the magnetic resistivity, and the unknowns are the fluid 
velocity field u(x, = (ui,U2,M3), the fluid pressure p(x, t), the magnetic field 
/B(x, = (61,82,63), and the magnetic pressure g(x, i), where x = {xi,X2,X3), 
and < > 0. Note that, a posteriori, one can derive that Vq = 0. Due to the 
fact that these equations contain the three-dimensional Navier-Stokes equations for 
incompressible flows as a special case (namely, when B = 0), the mathematical 
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1. The Inviscid Resistive MHD-Voigt Model 



The magnetohydrodynamic equations (MHD) arc given by 
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theory is far from complete. For a derivation and physical discussion of the MHD 
equations, see, e.g., [12]. For an overview of the classical and recent mathematical 
results pertaining to the MHD equations, see, e.g., [18, 20]. In this paper we study 
the inviscid = 0) case, with the following inviscid regularization of (1.1), 



(1.2a) -a^c'tAu + dtu + (u • V)u + V{p + -\B\^) = (B ■ \/)B, 

(1.2b) dtB+{u-V)B-{B-\7)u + \7q = nAB, 

(1.2c) V-B = V-u = 0, 

(1.2d) (u,6)|t=o = (uo,6o), 



where a > is a regularization parameter having units of length, and fi > 0. 
Note that when a = 0, we formally retrieve the inviscid {v = 0) resistive system 
(1.1). Furthermore, adding a forcing term to (1.1a) or (1.1b), or reintroducing a 
viscous term vAu to the right-hand side of (1.2a) (with u > 0), does not pose 
any additional mathematical difficulties to the results or the analysis of the present 
work, so it will not be discussed further. 

In [35], we studied inviscid, irresistive {v = 0, ^ ^ 0) magnetohydrodynamic 
equations with an additional inviscid regularization on both the velocity and the 
magnetic terms. That is, we studied the system 

(1.3a) -a^dtAu + dtu + (u • V)u + V{p + ^\B\^) = [B ■ V)B, 

(1.3b) -aljdtAB + dtB + (u • V)i3 - {B ■ V)u + = 0, 

(1.3c) V • S = V • u = 0, 

where a,aM > 0, the boundary conditions were taken to be periodic, and we also 
required 

(1.4) [ udx = I Bdx^O. 

Jn Jn 

Following the ideas of treating the Voigt- regularization of the 3D Euler equations, 
i.e., the inviscid simplified Bardina model, presented in [8], we proved in [35] that if 
Uo, i3o G H^{fl) and are divergence free, then (1.3) has a unique solution (u, B) lying 
in C^((— cx), cx)), 7J^(r2)). It should be noted that a related regularization known 
as the simplified Bardina model for the MHD equations has been studied in [33] in 
the viscous, resistive (j/ > 0, ^ > 0) case, with periodic boundary conditions. This 
model coincides formally with (1.3) in the inviscid, irresistive (:/ = 0, /i = 0) case, 
a phenomenon which was first explored in [8] in the context of the Euler equations. 
We also note that the viscous, irresistive {v > 0, /.t = 0) Bardina model for the 
MHD equations in was studied in [10]. 

The simplified Bardina model first arose in the context of turbulence models 
for the Navier-Stokes equations in [36]. Based on this work, the authors of [8] 
studied the simplified Bardina model in the special case where the filtering is given 
by inverting the Helmholtz operator / — a^A. In [8], the authors also studied the 
inviscid simplified Bardina model with this choice of filtering, which is known as 
the Euler- Voigt (or sometimes Euler- Voight) model, and proved the in their paper 
global regularity of solutions to the Euler- Voigt (i.e., inviscid simplified Bardina) 
system. Higher-order regularity, including spatial analytic regularity, of solutions 
to the Euler- Voigt model was then established in [35]. 
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The particular regularization considered in the present work, known as a Voigt- 
regularization, belongs to the class of models known as the a-models, which have a 
rich recent history (see, e.g., [8, 13-16, 24, 27, 28], also see, e.g., [8, 35] for historical 
discussions). Voigt-type regularizations in the context of various hydrodynamic 
models have been the focus of much recent research, see, e.g., [8, 10, 11, 21, 29- 
31, 31, 34, 35, 37, 46]. See also [21] for the application of Navier-Stokes-Voigt model 
in image inpainting. Voigt-regularizations of parabolic equations are a special case 
of pseudoparabolic equations, that is, equations of the form Mut + Nu = /, where 
M and N are operators which could be non-linear, or even non-local. For more 
about pseudoparabolic equations, see, e.g., [4, 9, 19, 45, 52-56]. 

We note that if one reintroduces a viscous term vAvl to the right-hand side of 
(1.2a) (with J/ > 0), it is possible to make sense of Dirichlet (no-slip) boundary 
conditions u = 0, and all of the results of the present work hold for such a case. 
Observe that the results reported here are also valid in the whole space R^^ by 
employing the relevant analogue tools for treating the Navier-Stokes in the whole 
space. We also remark that in the case > 0, /i > 0, and a > 0, one can consider 
the case of physical boundary conditions for (1.2), that is, uj^^ = 0, n • = 
and n X (V X ;B) Ian = 0, which are also known as no-slip, superconductor boundary 
conditions (see, e.g., [5, 20, 50]). With these boundary conditions, one can also 
prove that the system (1.2) enjoys global regularity with appropriate modifications 
to the methods employed in the present work, a subject of a forthcoming paper. 
However, the higher-order regularity of solutions to (1.2) in the case of physical 
boundary conditions does not follow directly from the proofs below, and one would 
have to modify the functional spaces taking into consideration the presence of a 
physical boundary (see, e.g., [32] and references therein). 

System (1.2) was introduced and studied in the two-dimensional case in [43], 
where global well-posedncss was proven under the assumption that Uq & V , Bq € 
H. The three-dimensional case was studied in [11], where global regularity was 
proven, assuming the initial data Uq S H^{^1), Bq G H^{n) and are divergence 
free. Here, we relax the hypotheses of the theorems given in [11] by requiring only 
that Uq € H^{n), Bq G L^(il) and arc divergence free for the existence of weak 
solutions, and uo,So G i/^(0) for the existence of strong solutions. Furthermore, 
we prove the uniqueness of strong solutions to (1.2), a result which is stated, but 
not proven, in [11]. We also prove the higher-order regularity of (1.2). In Section 
2, we introduce some notation and preliminary results. In Section 3, we establish 
the global existence of weak and strong solutions, and the uniqueness of strong 
solutions among the class of weak solutions. We give a complete proof using the 
Galerkin method, and we justify rigorously the a priori estimates and the existence 
and uniqueness results. In Section 4, we establish higher-order regularity of strong 
solutions. In Section 5, we show that strong solutions to the Voigt-regularized MHD 
equations (1.2) converge, as a — > 0, to strong solutions of the MHD equations (1.1) 
(with fj, > 0) on the time interval of existence of the latter. Furthermore, in 
Section 5, we establish a blow-up criterion for solutions of the inviscid, irresistive 
MHD equations, which can be easily implemented in numerical simulations. 

2. Preliminaries 

In this section, we introduce some preliminary material and notations which are 
commonly used in the mathematical study of fluids, in particular in the study of 
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the Navier-Stokes equations (NSE). For a more detailed discussion of these topics, 
we refer to [17, 25, 57, 58]. 

In this paper, we consider only periodic boundary conditions. Our space of test 
functions is defined to be 

V := 1^ e J" : V • ^ = and J ip{x) dx = o| , 

where T is the set of all three-dimensional vector-valued trigonometric polynomials 
with periodic domain = := [0,1]^. We denote by and iJ™ the usual 
Lebesgue and Sobolev spaces over il, and define H and V to be the closures of V 
in and H^, respectively. We define the inner products on H and V respectively 

by 

(u,v)='V / lijVjdx and ((u,v))= V / — ^— ^dx, 

and the associated norms |u| = (u, u)^/^, ||u|| = ((u, u))^/^. Note that ((•, •)) is a 
norm due to the Poincare inequality, (2.9) below. We denote by V the dual space 
of V . The action of V on V is denoted by (•, •) = (•, ■)y,. 

Let X be a Banach space with dual space X' . We denote by iP((a, 
the space of Bochner measurable functions t i~> w(i), where w(f) G X for a.e. 
t G (a, 6), such that the integral J^'' ||w(t)||^ dt is finite (see, e.g., [1]). A similar 
convention is used for C^((a^h\X). The space Ciu([0,r],X) is the subspace of 
i°°([0,r],X) consisting of all functions which are weakly continuous, that is, all 
functions w G L°°([0,r],X) such that (vif(i),v) is a continuous function for all 
V G X' . Note that here and below, we abuse notation slightly, writing w(-) for the 
map t H- vv(t). In the same vein, we often write the vector-valued function w(-,t) 
as w(t) when w is a function of both x and t. We stress that whenever we write 
Vif satisfies ^ G i^((a, we implicity mean that w : (a, 6) — !> X is absolutely 

continuous in time with values in X (see, e.g., [58]). 

Wc denote hy Pa ^ H the Leray-Helmholtz projection operator and define 
the Stokes operator A ;= —P^A with domain 'D{A) := {^V . A can be extended 
as a bounded linear operator A -.V ^ V, such that ||^v||y' = ||v|| for all v G T^. It 
is known that A^^ : H T^{A.) is a positive-definite, self-adjoint, compact operator, 
and that there is an orthonormal basis {wi},^-^ of H consisting of eigenvectors of 
A corresponding to eigenvalues such that Av/j = XjWj and < Ai < 

^2 < A3 < • ■ • (see, e.g., [17, 57, 58]) repeated according to their multiplicity. Let 
Hm ■= span {wi, . . . , w„j}, and let P,„ : H — >■ Hm be the orthogonal projection 
onto Hm with respect to Notice that in the case of periodic boundary 

conditions, i.e., in the torus T'^, we have A = —A, and Ai = {2tt)~'^ (see, e.g., 
[17, 57]). We have the continuous embeddings 

(2.1) V{A) '^V ^ H = H' ^V'. 

Moreover, by the Rellich-Kondrachov Compactness Theorem (see, e.g., [1, 22]), 
these embeddings are compact. 

It will be convenient to suppress the pressure term by applying the Leray- 
Helmholtz projection P^ and use the standard notation for the non-linearity, 

(2.2) S(u,v) :=P,((u.V)v) 
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for u, V G V. We list several important properties of B which can be found for 
example in [17, 25, 57, 58]. The proof of this lemma relics mainly on (2.2) and 
inequalities of the type (2.6)-(2.9) below. 

Lemma 2.1. The operator B defined in (2.2) is a bilinear form which can be ex- 
tended as a continuous map B : V xV ^ V' . Furthermore, the following properties 
hold. 

(i) For u, V, w € V , 
(2.3) (S(u, v), w)^, = - (B(u, w), v)^, , and (B(u, v), v)^,, ^ 0. 
(ii) We have the following estimates. 



(2.4a) 


(i3(u. 


v) 




<C|ur/2||u||V2||v|I||w||, 


V u e l^, V G w e y. 


(2.4b) 


(S(u, 


v) 




<C||u||||v|||w|V2||w||V2, 


V u e y, V e w e y, 


(2.4c) 


(S(u 


v) 




<C|u|||v||V2|Av|V2||w||, 


Vuei/,vGX'(A),wey, 


(2.4d) 


(B(u 


v) 




<C||u||||v||V2|Av|i/2|w|, 


Vu e y,v e P(A),w e i?. 


(2.4e) 


(B(u 


v) 




<C||ur/2|Au|V2||v|||w|, 


Vu e V{A),^^ e V,w eH, 


(2.4f) 













<C|u|V^||ur/^|v||Aw|, 



V2||u||V2| 

< C|uPv||w|i/2||w|li/2, 



Vuey,veiJ,we P(^), 
VueiJ,veX>(A),wGy. 
by 



I (5(u,v),w)p(^), 
(2.4g) |(S(u,v),w)^, 

We also define the trilinear form b : V x V x V 

(2.5) 6(u,v,w) := (B(u,v),w)^,. 

Here and below, Ka, Ka,^, C{- ■ ■), etc. denote generic constants which depend only 
upon the indicated parameters, and which may change from line to line. Let us 
distinguish between K and C. K will depend on some norm of the solutions, but 
C will not be used for constants which depend on functions. 

Next, we recall Agmon's inequality (see, e.g., [2, 17]). For u G 'F>{A) wc have 

(2.6) ||u|U^(a)<C||u||V2|Au|i/2 . 

We also have the Sobolev and Ladyzhenskaya inequalities in three-dimensions, 

(2.7) ||u|U3<qu|i/2||u||V2 

(2.8) ||u|U. <C||u||, 

for every u £ V. Furthermore, for all w G y, wc have the Poincare inequality 

(2.9) |w| < A^'/'|Vw| = A7'/'||w||. 

Due to (2.9) and the elliptic regularity of the solutions to the Stokes equation (see, 
e.g., [17, 57]), for w G ^{A), we have the norm equivalence 

(2.10) |Aw| ^ ||w||^2. 

Finally, we note a result of deRham (see, e.g., [58, 59]), which states that if g is a 
locally integrable function (or more generally, a distribution), we have 

(2.11) g = Vp for some distribution p iff (g, v) =0 Vv G V. 

This result is useful for recovering the pressure term as it is treated, for example. 
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3. Existence and Uniqueness of Solutions 

This section is devoted to stating and proving our main result. As mentioned 
in the introduction, we relax the conditions of the results of [11], where it was 
assumed that Uq G V{A), Bq € V to derive the existence of a strong solution. 
Here we define the notion of a weak solution to (1.2), for which we only need to 
assume that Uq ^ V , Bq ^ H to show global existence without uniqueness. As for 
a strong solution, we only need to assume Uq,Bo £ V to prove global existence, 
uniqueness, and continuous dependence on initial data. We note that, although the 
major features of this model that allow for a proof of well-posedness were exploited 
formally in [11], the a priori estimates can be sharpened. This allows us to define 
the notion of weak solutions and to prove their global existence. We also prove the 
uniqueness of strong solutions, which was stated without proof in [11]. Furthermore, 
there are subtleties in passing to the limit due to the addition of the Voigt term 
— a^(?t Au in the momentum equation (1.2a) that one has to address in the rigorous 
proof. Here, we give a fully rigorous derivation and justification of these estimates, 
as well as the passage to the limit. 

In order to prove global existence, we use the Galerkin approximation procedure, 
based on the eigenfunctions of the Stokes operator. The proof is broken into several 
steps. First, we show that there exists a solution to the finite-dimensional approxi- 
mating Galerkin problem which is bounded in the appropriate norms. In particular, 
we show that the time derivatives of the sequence of approximating solutions are 
uniformly bounded in the appropriate spaces. We then extract appropriate sub- 
sequences using the Banach-Alaoglu and Aubin Compactness Theorems (see, e.g., 
[17], p. 68-71 or [47, 58]), and pass to the limit to obtain a global solution to system 
(1.2). Finally, we argue that the solutions satisfy the initial conditions in the sense 
given in Definition 3.1 below. 

3.1. Existence of Weak Solutions. Before we begin, we rewrite (1.2) in func- 
tional form. Applying Pa to (1.2) and using the notation introduced in Section 
2, we obtain the following system, which is equivalent to (3.1) (see, e.g., [58] for 
showing the equivalence in the context of the Navier-Stokes equations) 

(3.1a) ^ (a^Au + u) = B(B,B) - B(u,u), 

dt 

(3.1b) ^B + ^lAB ^ B(B,u) - B(u,B), 

dt 

(3.1c) B{Q) = Bo, u(0) = uo, 

where (3.1a) is satisfied in the sense of L^/^{{0,T),V'), (3.1b) is satisfied in the 
sense of L'^{{0,T),V'), and (3.1c) is satisfied in the sense of Definition 3.1 below. 
Systems (1.2) and (3.1) are equivalent, and one can recover the pressure terms p and 
q {q = 0) by using (2.11), as it is done for the case of the Navier-Stokes equations 
(see, e.g., [20, 58]). 

Definition 3.1. Let Uq £ V, Bq £ H. We say that (u, B) is a weak solution to (3.1), 
on the time interval [0,T], if u e C([0, T], F), B € L^{{0,T),V) D Cr,{[0,T], H), 
^ G L^{{0,T),H), ^ e L^{{0,T),V'), and furthermore, (u,S) satisfies (3.1a) in 
the sense of ^^/^([O, T], V) and (3.1b) in the sense of ^^([o, T],V'). Furthermore, 
if Uo, Bo G y, we say that (u, B) is a strong solution to (3.1) if it is a weak solution. 
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and additionally, B G L^{{0,T),V{A)) D C{[0,T],V), ^ G C{[0,T],V), and ^ 6 
i2((0,r),i7). 

With this definition, we are now ready to state and prove the following theorem. 

Theorem 3.2. Let Uq E V, Bq ^ H. Then (3.1) has a weak solution (u, ;B) for 
arbitrary T > 0. 

Proof. Let T > be fixed. Consider the finite dimensional Galerkin approximation 
of (3.1), based on the eigenfunctions of the operator A (see Section 2), given by the 
following system of ODEs in H„i x Hm- 

(3.2a) — (um + a^Aum) + P,„_B(u )=P,nB{B,n,B,n), 

(3.2b) ^B,n + tiABra + P™B(u„„ Bra) = PmB{Bra. U™), 

at 

(3.2c) Z?,„(0) = P,„So, u,„(0) = P™uo. 

We look for a solution u^, Bm G C-^([0, T^), i?m) of (3.2). By applying the operator 
(/ + a^A)^"^ to (3.2a), we see that (3.2) is equivalent to a system of the form 
y = F(y), where F : H„i x Hm Hm x Hm is a quadratic polynomial. By classical 
ODE theory, this system has a unique solution on [0,Tm) for some > 0. Let 
[0,T^'^^) be the maximal interval where existence and uniqueness of the solution 
of (3.2) holds. 

Next, we show that T^'^^ = oo. Indeed, taking the inner product of (3.2a) 
with Um{t) and (3.2b) with Bm{t), for t £ [0,T^^^), and integrating by parts with 
respect to the spatial variable and using (2.3), we have. 



(3.3a) ~ («'||u,„f + |u,„|2) = (B(B„„S,„),u„), 

(3.3b) il|S,„|2+^||^,„||2 ^ {B{Bm,Urn),Bm) ^ -{B{Brn,Bm),Urn). 

Adding (3.3a) and (3.3b) gives 

(3.4) i| (a'l|u™|p + |u,„p + |B„,p) = -Ml|S™f < 0. 



Integrating the equality in (3.4) in time, we obtain for t G [0,T'^^^) 

a2||u„(i)f + |u„(<)|2 + |Z?,„(t)|^+2Ai / \\B,n{s)fds 

Jo 

= a2||u„(0)f + |u,„(0)p + |B„,(0)p 
(3.5) < (Kif :=a2||uof + |uop + |So|'. 

This bound, together with the fact that the vector field F(y) in (3.2) is a quadratic 
polynomial, imply that T™^^ = oo. Furthermore, we see from (3.5) that for fixed 
but arbitrary T > 0, we have 

(3.6a) u„ is bounded in L°°([0,r], V), 

(3.6b) Bm is bounded in L°°{[0,T],H) n L\[0,T],V), 

uniformly with respect to m. 

As mentioned at the beginning of this section, our goal is to extract subsequences 
of {u„i} and {Bm} which converge in -L^((0, T),H) by using the Aubin Compactness 
Theorem (see, e.g., [17], p. 68-71 or [47, 58]). To satisfy the hypotheses of Aubin's 
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theorem, we show that ^ is uniformly bounded in L^{{0, T), H) ^ L^{{0, T),V'), 
and that is uniformly bounded in L^((0, T), 1/'), with respect to m. Using 
equation (3.1a), we have from (2.4f), 



dt 



ViA)' 



< \\P,nB{B,n, B,n)\\v(Ay + || ^'mB(u,„ , U„i) || •p(^)/ 

(3.7) < CiKif/^WBrnW'/'' + {Ki)'a-'/^ 

where we have used (3.5). Estimating differently, we have 



(/ + a^A) 



dUr, 



dt 



< \\PmB{B,n,B„,)\\v' + \\P,nB{u^,Urn)\\v' 

= sup {B{B,n,B„i),P„iw) + sup (B(u„,u„,),F„iw) 

||w|| = l |[w|| = l 

<C sup |S„,|l/2|16„f/2||w||+C sup |u,„|l/2||u„||3/2|l^|| 
||w|| = l ||w|t = l 

(3.8) <C{KinB„M'/' + {Kifa-'/^ 

Thus, due to (3.5), the right-hand side of (3.7) is uniformly bounded in L^(0, T), and 
right-hand side of (3.8) is uniformly bounded in L''/^(0, T), and hence, (I+a'^A)^^ 
is uniformly bounded in L''([0, T], 23(A)') n L'^/^{[0,T],V') with respect to 771. By 
inverting the operator (/ + a^A), we have 

(3.9) 



dt 



is bounded in L*([0,r],iJ) and L'^/^{[0,T],V), 



Next, we estimate From equation (3.1b) we have, thanks to (2.4a) and 



uniformly with respect to m. 

Nex 
(2.4b), 

-7^ < mnB{B,n,U.^)\\v' + \\Pr,rB{u^,B^)\\v, + ^l\\ABm\\v' 

dt y, 

<C\B„,\'/^\\B„^\\'/^\M\+^i\\Br,^\ 

(3.10) <C(A-i)3/2a-i||S„||V2 + ^||S„||, 

where the last estimate is due to (3.5). Thus, by virtue of (3.6b), it follows that 



(3.11) 



dB„. 
dt 



is bounded in L^{[0,T],V'), 



uniformly with respect to m. 

We have shown in (3.6a) that u„i is uniformly bounded in L°°{[0,T],V) ^ 
L2([o,T],F), that ^ is uniformly bounded in L^{[0,T],H) ^ L^{[0,T],V'), that 
i3™ is uniformly bounded in L°° {[0,T], H)nL'^{[0,T],V), and that ^ is uniformly 
bounded in L'^{[0,T],V'). Thus, by the Aubin Compactness Theorem (see, e.g., 
[17], p. 68-71 or [47, 58]), there exists a subsequence of {Bm,Um) (which we relabel 
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as (i3m,Um), if necessary) and elements B,uE L'^{[0,T],H) such that 
(3.12a) Bm^B strongly in L2([0,T],i7), 

(3.12b) Um^u strongly in L2([0,T],i7). 

Furthermore, using (3.6a), (3.6b), (3.9), (3.11) and the Banach-Alaoglu Theo- 
rem, we can pass to additional subsequences if necessary (which we again rela- 
bel as {Bm,nm)), to show that, in fact, u € L°°([0, T], F), B G L°°([0, T], 7J) n 
L%[0,T],V), iuGL^i[0,T],V)n i4/3 ( [0, T] , if ) , As 6 L2 ( [0, T] , r ) , and 

(3.13a) Bm ^ B and ^ u weakly in L^([0, T], F), 
(3.13b) Bra^B and u„ ^ u weak-* in i°°([0, T], iJ), 
(3.13c) u„^u weak-* in i°°([0,r],F), 

(3.13d) -^u„ ^ -^u weak-* in L''([0,T],i?) and L''/3Qo,T],F), 

(3.13e) ^Bm^Y^ weak-* in L2([o,T],F')- 

as m ^ oo. 

Let k be fixed, and take m > k. Let w e C'^{[0,T], Hk) with w(T) = be 
arbitrarily given. By taking the inner product of (3.2) with w, integrating on 
[0,r], and using integration by parts, we have 

(3.14a) - (u,„(0),w(0)) - a2((u,„(0),w(0))) 

{u,n{t),w'{t))dt + a^ / ((u,„(t),w'(t)))di 



{B{Bm{t),Bra{t),Pra^{t))dt- [ (S(u,„ (t) , U„ (t)) , P,„ w(t )) , 



(3.14b) -(B„,(0),w(0))- / iBUt),Mt))dt + ii I ((6„,(t),w(t)))dt 



iBiB,nit),U„,{t)),P„Mt))dt- / iB{u,nit),B„,{t)),P,r^w{t))dt. 

Jo 

We show each of the terms in (3.14) converges to the appropriate limit, namely, we 
will find that equations (3.14) hold with {Bm, Um, Pm} replaced by {B, u, /}, where 
/ is the identity operator. First, thanks to (3.13a), we have 

r [[B^{t)Mt)))dt^ 11 r {{B{t)Mt)))dt, 
Jo Ja 

(u„(t),w'(t))dt ^ / {u{t),w'{t))dt, 



a" ((u™(t),w'(t)))dt ^aM ((u(t),w'(t)))dt. 



T 



{B„,{t),W(t))dt^ / {B(t),W(t))dt 



Next, we must show the convergence of the trilinear forms. We will only show 
the convergence of one of them, as the rest are similar (see, e.g., [17, 58] for similar 
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arguments in the case of the Navier-Stokes equations). Namely, we wiU show that 

/(m):= / {B{u„,{t),u„,{t)),P„Mt))dt~ f {B{u{t),u{t)),w{t)))y, dt ^ 0. 
Jo Jo 

To this end, let 

him) := [ {B{ur,^it) - u(t),P™w(t)),u™(t))^, dt, 
Jo 

him) := [ (B(u(i),P,„w(t)),u,„(t) -u(i))^, dt, 
Jo 

and note that /(m) = Iiim) + him), where we have used (2.3). Since w £ 
C^([0,T],fffc) and fc > to, we have P,„w = w. Thus, thanks to (2.4c), (2.4d), and 
(3.6a), a simple application of Holder's inequality and (3.12b) shows that him) 
and him) — > 0, and thus /(to) -> for w G C^([0, T], iJ^). As mentioned above, 
similar arguments hold for the other tri-linear terms. Note that Um(0) := -PmUo — > 
Uq in V and ;B„i(0) := PmSo — > in H. Thus, passing to the limit as to — ?> oo in 

(3.14) , we have for aU w e ([0, T], fffc) with w(T) = 0, 

(3.15a) - (uo, w(0)) - «2((uo, w(0))) 

iu{t),w'{t))dt + a' f ((u(t),w'(t)))dt 
Jo 

{BiBit),B{t),wit))dt- [ (B(u(t),u(t)),w(t))di, 
Jo 

(3.15b) -{Bo^wiO))- [ {B{t),w'it))dt + ^, [ ((B(t),w(0))di 

Jo Jo 

'^{BiBit), uit)), w(0) dt - riB{uit), B{t)),wit)) dt. 
10 Jo 

Since Ci([0, T], iJ^) is dense in C^{[Q,T],V), we use (2.4a) and the facts that u G 
L°°((0,T),\/) and B £ L°°((0, T), i/) n ^^((O, T), y) to show that (3.15) holds for 
all w e C^i[Q,T],V) with w(T) = 0. In particular, (3.1a) and (3.1b) are satisfied 
by (u, S) in the sense of V , where the time derivatives are taken in the sense of 
distributions on (0,T). Allowing (3.1a) and (3.1b) to act on w and comparing with 

(3.15) , one finds that u(0) + a^AuiQ) = Uq + o^Auq and S(0) = So (see, e.g., [58, 
p. 195]). By inverting / + a^A, wc then have u(0) = uq. 

We must now show that u and B satisfy the requirements for continuity in time 
in Definition (3.1). Taking the action of (3.15b) with an arbitrary v S V and 
integrating in time, we obtain, for a.e. U),ti e [0,T], 

(3.16) (S(ti) - S(to), v) + r UBit), v)) dt 

Jto 

ti rti 



(i?(B(t),u(t)),v)^, dt- / {Biu{t),Bit)),w)y, dt. 

Jto 

Since the integrands are in Li((0, T)), we see from (3.16) that B G Cu,([0, T], V) by 
sending ti to. By the density of V in H and the fact that B £ L°°([0, T], iJ), a 
simple application of the triangle inequality shows that B G Cui{[0,T], H). Next, 
by (3.9), we have f^u € L'^i[0,T],H) ^ L^i[0,T],V'). Since we also have u e 
L2([0,T],F), it follows that u G C([0,r],iJ) (sec, e.g., [47, Corollary 7.3]). We 
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have now shown that u and B satisfy all the requisite conditions of Definition 3.1, 
and therefore we have proven the global existence of weak solutions. □ 

3.2. Existence of Strong Solutions. Next, we prove the existence of strong 
solutions, assuming Bq G V. Uniqueness of strong solutions among the class of 
weak solutions, and also higher-order regularity will be proven in the next two 
sections. 

Theorem 3.3. Let Uq,Bo G V. Then for every T > 0, (3.1) has a strong solution 
(u,i3) on [0,r]. 

Proof. Taking the inner product of (3.2b) with ABm, we use (2.4a) and (2.4e) to 
obtain, 

il\\B^\\^ + ^\AB.^\^ 

= -{B{Bm,Urn),AB„,) + {B{u,n,B^),AB,n) 
<C\\Brn\\'/MBrn\'^'\\Urr,]\\ABrn\+C\\UrJ\\^^^^^^ 

(3.17) < 

since ||u„i|| is uniformly bounded by (3.5). Due to Young's inequality, we have 

(3.18) a-iifi||i3„r/2|AS„|3/2 < KlJB„r + ^\AB^\\ 

where K'^ ^ := C(a^^i^„)^^~^. Therefore, combining (3.17) and (3.18), we have 

(3.19) + ^^1^^™!' ^ Klj^rr^f. 

Integrating (3.19) on [0,t] gives 

\\B,^{t)f+iil \ABm{s)\^ds<\\BUOW + ^Kl^ f ||B„(s)fds 
^0 Jo 



(3.20) <\\Bo\\^ + Kl^^:=Kl^ 



due to (3.5). Since we are now assuming Bq G V (3.20) implies 
(3.21) B^ is bounded in L°°([0,T],F) ni2([0,r],2?(A)). 

uniformly with respect to m. Furthermore, recalling (3.7), the improved bound 

(3.21) now yields 

(3.22) ^ is bounded in ( [0, T] , F) . 

at 

uniformly with respect to m. 
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Next, we estimate From (2.4a), (2.6), (2.9), (3.1b), and (3.5), we have. 



dB„ 



dt 



dt 

< \B{B,n,u„,)\ + |B(u,„,i3,„)| +Ai|^S,„| 

= sup (B(B™,U™),w) + sup {B{u.,n,B,n),w) + ^\ABr,i\ 
|w| — 1 |w| — 1 

< sup ||B™||l~|Vu,„||w| + sup |u,„|i/2||u„l|i/2|^^^||^|+^|^^^^| 

|w| — 1 jw| — 1 

< C|iS„jr/2|AS„Ji/2||u„|| +/i|Ai3™| 
Thanks to this and (3.6b) 

(3.23) ^ is bounded in ^^([o, T], iJ), 

dt 

uniformly with respect to m. Now, from the proof of Theorem 3.2, we aheady know 
that there exists a weak solution (u, B) of (3.1) such that Bm' ;B in i°°([0, T],H) 
and L'^{[0,T],V) for some subsequence {Bm'}- Thanks to (3.21) and (3.23), we 
may apply the Aubin Compactness Theorem (see, e.g., [17], p. 68-71 or [47, 58]) to 
extract a subsequence (relabeled as (u„i,_Bm)) such that 

(3.24) B,n B strongly in ^^([o, T],V). 

Using (3.21), the Banach-Alaoglu Theorem, and the uniqueness of limits, we 
may pass to additional subsequences if necessary to show that B G L°°{{0, T),V)r\ 
L2((0,r),2?(A)). It is easy to see from (3.23) that ^ € L^{{0,T),H). Thus, we 
must have B G C([0,T],F) (see, e.g., [47, Corollary 7.3]). Finally, since u,B e 
C{[0,T],V), it follows easily that the right-hand side of (3.1a) is in C{[0,T],V'). 
Inverting {I + a^A) shows that ^ € C{[0,T],V). Therefore, we have shown the 
existence of a strong solution to (3.1). □ 

3.3. Uniqueness of Strong Solutions and Their Continuous Dependence 
On Initial Data. In this section, we prove the uniqueness of strong solutions 
among the class of weak solutions. As mentioned in the introduction, the uniqueness 
of (strong) solutions is stated in [11], but no proof is given. We begin with a lemma, 
which is reminiscent of the Lions-Magenes Lemma (see, e.g., [39, 58]). We use 
similar ideas to those in [58, Lemma 1.2, p. 176]. 

Lemma 3.4. Let v e C((0,T),7J) and S LP{{Q,T),H) for some p g [l,oo]. 
Then the following inequality holds in the in distribution sense on (0,r). 

Moreover, ]v]^ is absolutely continuous. 

Proof. First, we note that (3.25) makes sense, due to the fact that 1 1— J- lv(<)p and 
t I— >■ (^v(t),v(<)) are both elements of L^{[0,T]). Let us write v for the function 
which is equal to v on [0, T] and equal to on R\ [0, T]. By a standard mollification 
process, we can find a sequence of functions {v/i;}j.gj^ in C°°([0, T]) such that v 
in Lj^^((0,r),i7) and ^v,. ^ in Lf^ ^((0, T), 7J). Clearly, equality (3.25) holds 
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for Vfc, and also |vfcp |v|2 in ^^^^((O, T)). We also have for <ti <t2 <T, 
I (^v., V,) - ( Av, v) I < r K^v, - ^v, V,) I + I ( A V, V, - v) I 
^ '^Wli^k - ft^\\Lpat,.t2))\Wk\\L--iit,,t2)) + \\ft^\\LpHt,,t2))\Wk - v|lLoo((t^,t^)). 

Thus, (^v/j,v/j;) — > (^v,v) in Lj'Q^((0, T)) as wcU, and so (3.25) holds in the scalar 
distribution sense. Furthermore, since the right-hand side of (3.25) is integrable, 
we have |vp G W^'^{{0,T)), and so |vp is absolutely continuous in time. □ 

Theorem 3.5 (Uniqueness and Continuous Dependence On Initial Data). Let 
{u^,B^) be a strong solution to (1.2) with initial data Uq, Bq € V and let (u^,B^) 
be a weak solution with initial data Ug, Bg G V . Let us write <5v := — for two 
arbitrary, consecutively labeled vectors and . We have, 

(3.26) |5u(i)|'+a'||<5u(t)f + |5e(i)|2+;, / \\6B(s)fe'^'-~^^^'-'Us 

Jo 

< {\duo\^ + a^\\Suo\\^ + \6Bo\^) e<-* 

= fi{K\ ^, ^,a,fi) is a positive constant. In particular, if Uq = Uq and 
Bo s Bq in the sense of H, then = and B^ = B"^ in the sense of L^{[0,T],H). 

Proof. Let us denote Su := — and 6B := B^ — B^. Subtracting, we get an 
equation for 5u: 

(3.27) (/ + a'^)^ = i^S ■ ^)B^ + (S' • ^)5B - (d^u • V)u^ - [v? ■ V)5n. 
Applying (/ + q;2^)-i/2 to both sides of (3.27) we obtain 

(3.28) [I + a^Af/^'^ 

at 

= {I + a^AY^I"^ {{5B ■ V)B^ + [B^ ■ V)SB ~ (Su ■ V)ui - (ua ■ V)^u) . 

Since <5u G C{[0,T],V) and ^ G L*/^{[0,T],V), and therefore (I + a^A)^/^^ G 
L''/3([0,T],7J), we may take the inner product of (3.28) with (/ + a^A)'^/'^6u G 
C([0, T],H). Thanks to Lemma 3.4, we have ((/ + a'^Ay/'^f^Su, {I + a^A^/'^Su) = 

1^ |(/ + Q,2A)i/25u| . After using (2.3), we arrive at 

(3.29) |(/ + a2A)i/25up = b{6B,B\Su)+b{B'\SB,5u)-b{6u,u\Su). 

Arguing in a similar way to the standard Navier-Stokes theory (in particular, using 
the Lions-Magenes Lemma [39, 58] rather than Lemma 3.4, see, e.g., [17, 47, 58]), 
we derive the following equation for SB. 

(3.30) + Mll^^f = K^B^ u\SB) + 6(S^ Su, SB) - 5((5u, B\SB). 

Since b{B^,Su,SB) = -b{B'^,SB,5u) by (2.3), we may add equations (3.29) and 
(3.30) to obtain 

i^(|(/ + aW/^<5u|' + |<5B|^)+Ml|5S||^ 

= b{SB, B\Su) - b{Su, ui , Su) + b{SB, u\SB) - b{Su, , SB) 
<C||<5S||||Si||||Ju||+C|H'/'ll'5uf/^||ui|| 
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+ C|(5S|i/2||5Bf /2||ui|| + C\\Su\\\\B^\\\SB\'^^\\dB\\^/^ 
<KlJSB\\\\Su\\+KiJSunSur/' 

+ <^|5B|l/2||5Bf/2+i^3 J|5u|||JB|l/2||5S||l/2 

<lKl,{\Su\^- + a^Sur + \6B\') + ll\\6Br, 

where K^^^ = K^^^{K^^^, K^^^,a, fi) is a positive constant. Here, we have used 
Young's inequahty several times, as well as (2.4a), (2.4b), (2.9), (3.5), and (3.20) 
(which hold in the limit as m — > cx) for ,B*,u' by properties of weak convergence). 
Therefore, 



dt 



2 



(/ + a'Af/^Sn + \5B\'' + ii\\5Bf < K^^ {\5uf + a^\\Suf + \SBf 



Using Gronwall's inequality together with the identity |(/ + a^A)^/^5u|'' — |(Sup + 
a^||(5u|p now implies (3.26). □ 

4. Higher-Order Regularity 

We now prove that the solutions to (3.1) (equivalently (1.2)) enjoy regularity 
whenever Uq, ;So G H'^ D V for s > 1. We note that the calculations in this section 
are done formally, but can be made rigorous by proving the results at the Galerkin 
approximation level, and then passing to the limit in a similar manner to the 
procedure carried out above (see, e.g., [35]). By the uniqueness of strong solutions, 
all strong solutions arc regular. 

Theorem 4.1. Consider (3.1) under periodic boundary conditions fl = T''. Sup- 
pose for s >l that uq e H'T] V, Bq e HT\ V. Then (1.2) (equivalently (3.1) ) has 
a unique strong solution (u,B) with u G L°°([0, T], if" DV), B € L°°{[0,T],H' n 
V) n L'^{[0,T], H''^^ nV). if we furthermore assume that uo G H'+^ n V, then we 
also have u G L°°([0, T], i/'^+i n F). 

Remark 4.2. It is possible to extend Theorem 4.1 to include the case of a specific 
type of Gevrey regularity, which is analytic in space, as was done for the Euler-Voigt 
equations in [35]. For the sake of brevity, we do not pursue such matters here. For 
more on Gcvrcy regularity, see, e.g., [6, 7, 23, 23, 25, 26, 29, 32, 38, 42, 44, 48] and 
the references therein. 

Proof. As indicated above, we only give a formal proof, but the details can be 
made rigorous by using the Galerkin approximation procedure (see, e.g., [35] for 
a complete discussion of this method applied to proving higher-order regularity in 
the context of the Euler-Voigt equations). Furthermore, we only prove the result 
in the cases s = 1,2, as the cases s > 2 are more complicated notationally, but 
not conceptually (see, e.g., [35]). In the case s = 1, the first statement has already 
been settled by Theorems 3.3 and 3.5. To prove the second statement, assume that 
Uo G 2?(^) and Bq G V. Let /3 be a multi- index with \(3\ — 1. Applying to 
(1.2a), integrating the result against d^u, and integrating by parts, we have 
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{a^\\d^u\\^ + |9^u|2) = {d^B-'7B,d^u) + {B ■ SId^B,df^M) 

- {{d^u ■ V)u), a^u) - {B{u ■ y)d'^u),d^u) 

(4.1) ={df^B-VB,d'^u) + {B-Vd^B,d^u) 

Where we have used (2.3). Next, we apply d'^ to (1.2b), integrate the result against 
d^B, and integrate by parts to find 

\^\d^l3\^ + Mll^'^^f = (d^S ■ Vu,d^B) + {{B ■ V)d^u, d^B) 

- ((a^u ■ V), S), a^S) - {{u ■ V)d^B), d^B) 

(4.2) = {d^B ■ Vu, d^B) + {{B ■ V)5'^u), dl^B) 

- {{df'u-V)B),d^B). 

Since {{B-V)d'^u,d^B) = -{{B-V)d^B, df^u) by (2.3), we may add (4.1) and (4.2), 
to obtain a very important cancellation of the terms that involve the highest order 
derivatives: 

1±{\0PB\' + a^dPur + \dPu\') + ^A\^PBr 
= {{d^B ■ V)B, 9'^u) - {{df^ ■ V)u, u, d^u) 

+ {{d^B ■ V)u, d^B) - {{d^ ■ V)u, B, d>^B) 
<C||S||||S||^3/.|Au|+C||u||||u||^:v.|Au| 

< C||ef /2|AS|i/2|^u| + C||uf /2|^u|i/2|Au| 

< C||Sf /2|A6|i/2|Au| + A'i|Au|3/2 

< C\\Bf/^\AB\^/^\Au\ + Xi|Au|3/2 

< ||AZ5p + C||Bf |Au|4/3 + Xi|Au|3/2, 

since ||u|| is uniformly bounded on [0,r]. Summing over all /3 with \/3\ = 1, we have 

(4.3) ^ + a2|Au|2 + ||uf ) + ^i|ABp < C\\Bf\I\xx\^l^ + A'^ |Au|3/2. 

Letting y = l + ||Bp+a^|Aup + ||u||^ and dropping the term /.t|A^?p for a moment, 
we have an equation of the form 

where K{t) = C ■ K^{1 + ||S(t)p). Gronwall's inequality gives 
\\B{t)r + a'\Auit)\^ + \\uit)r 

(4.4) < At \\Bof + a'\\Auof + ||uof + K {1 + \\B{s)f) d)j . 

Since \\B{s)\\^ ds < oo by Theorem (3.3), we see that u E L°°{[0,T],V{A)), 
thanks to (4.4) and the norm equivalence (2.10). This in turn implies that the 
right-hand side of (4.3) is finite on [0,r]. Integrating (4.3) on [0,T], we find that 
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B G L^([0, T], 2?(A)). Thus, we have formally established the theorem m the case 
s = 1. 

Let us now take s = 2. To begin, wc formally take the inner product of (1.2a) 
with A'^u (recalling that, in the periodic case, A = —A) and integrate by parts 
several times to obtain 

i-i(o'||A„f+|AuP) 

= (AS • VB, Au) + 2{VB ■ VVB, Au) + {B ■ VAB, Au) 

- (Au • Vu, Au) - 2(Vu • VVu, Au) - (u • VAu, Au) 

(4.5) = [AB ■ VB, Au) + 2{VB ■ VVB, Au) + {B ■ VAB, Au) 

- (Au • Vu, Au) - 2(Vu • VVu, Au). 

where we have used (2.3). Next, we take the inner product of (1.2a) with A'^B and 
integrate by parts several times to obtain 

= {AB ■ Vu, AB) + 2{VB ■ VVu, AB) + {B ■ VAu, AB) 

- (Au • VB, AB) - 2(Vu • VVS, AB) - (u ■ VAB, AS) 

(4.6) = (AS • Vu, AS) + 2(VS • VVu, AS) - (S • VAS, Au) 

- (Au • VS, AS) - 2(Vu • VVS, AS). 

where we have again used (2.3). 

To prove the first statement of the theorem with s = 2, we take uq. So G T^IA). 
We already have u G L°°{[Q,T],'V{A)) by the case s — 1, so it remains to show 
S G £°°([0, T],V{A)) n L2([0, T],H^r\ V). Estimating the right-hand side of (4.6), 
we find 

(4.7) ^^l^^l' + mII^^II' < (^(|AS|||AS|||Au| + 6|Au||AS|3/2j| AS||^/2) 

< (i^^)i/2(|AS|||AS|| + |ASp/2||Agj|i/2) 

<<,|ASp + |||ASf, 

where K^^ = K%^^{K%,^i). Here, we have used (2.4g), (2.4b), and (2.4c) in the 
first inequality, (4.4) for the second inequality, and Young's inequality for the last 
inequality. Rearranging and using Gronwall's inequality, we find, 

(4.8) |AS(t)p + /' ||AS(,s)|| V<-(*-^) ds < lASope^^"-*. 

"'0 

Thus, by eUiptic regularity, we have S G L°°{{0, T],V{A)) n ^^([o, T],H^n V). 

We now prove the second statement of the theorem with s = 2. This statement 
can be proven independently of the first statement, so we do not rely on (4.8) in 
the calculations below. Let uo E DV and So G 2?(A). Adding (4.5) and (4.6) 
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and again noting the important cancellation of higher-order derivatives, we obtain, 
i|(|ABp + a^||Au|p + |Aun+;.||AZ5f 

= (AB • Vu, AB) + 2{VB ■ VVu, AB) - (Au • VB, AB) - 2(Vu • VVB, AB) 
+ {AB ■ Vi3, Au) + 2(V6 • VVS, Au) - (Au • Vu, Au) - 2(Vu • VVu, Au) 

< C{-i\AB\^'^\\ABf'^\An\ + 3|Au|5/2|| AuH^/^ ^ Q\\B\\\Au\^''^\\Anf^\\AB\\) 

< C • (1 + /^^)2(|AS|3/2|[AS||i/2 + ||Auf/2 + IIAuf /2||AS||) 

< C • (1 + Klf{\AB\^'^\\AB\\''^ + II Auf /2 + || Auf as||) 

< • (1 + Klf{\AB\^ + II Au|| + II Au||) + ^11 ABf . 

The first inequality is due to (2.4). The second inequality is due to (4.4). The third 
inequality follows from the norm equivalence (2.10), and the last inequality is due 
to Young's inequality. This leads to 

^ (|AB|2 + c.2||Au||2 + |Au|2) + ;i|l ABf 

<<,,(1 + |AB|2+C.2||AU||2+|AU|2), 

where — Ca.i_i ■ (1 + Ka)^ ■ Employing Gronwall's inequality, we have 
|Ai3(i)|2+a2|[Au(t)f + |Au(t)|2+^ / ||ABfe^=.''(*-^' 

"'0 

< (1 + lABoP + a^llAuof + lAuoHe^"-* 

Due to the equivalence of norms given in (2.10), we see from the above inequality 
that B e L°°{{0,T),V{A))nL^{{0,T),H^r]V) (which was obtained independently 
in (4.8)), and u e L°°((0, T), n V). Thus, we have formally proven the theorem 
in the case s = 2. A similar argument can be carried out inductively for the cases 
s > 2. See, e.g., [35] for a rigorous, detailed discussion of this type of argument. □ 

Remark 4.3. Note that one may add a diffusion term z^Au to the right-hand side of 
(1.2a) or a suitable smooth forcing term to one or both equations (1.2a) and (1.2b) 
and recover the same existence, uniqueness, and regularity results obtained above 
using only slightly modified techniques. Furthermore, doing so allows one to study 
the attractor of such a system, which has been done in the case of the Navier- 
Stokes-Voigt equations in [29, 30]. We will study these ideas in a forthcoming 
paper. 

5. Convergence and a Criterion for Blow-up 

In this section, we prove that solutions of the regularized system (1.2) converge to 
solutions of the non-regularized system (that is, the equations obtained by formally 
setting a = in (1.2)). We first state a short-time existence and uniqueness theorem 
which was proven in a slightly more general context in [3] (see also [49, 51]). 

Theorem 5.1. Given initial data Uo,;Bo G H V, there exists a time T > 
and a unique solution {u,B) to the system (1.1), with v = 0, fi > 0, such that 
u,Be C{[0,T],H^ nv). 
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Theorem 5.1 can be proven similarly to the case of the 3D Euler equations (see, 
e.g., [40, 41]). For example the existence of solutions can be proven by considering 
the finite-dimensional Galerkin approximations to (1.1) based on the eigenfunctions 
of the stokes operator, and showing that the sequence of solutions as the dimension 
increases is a Cauchy sequence, and that it converges to a solution of (1.1) in 
an appropriate sense. Theorem 5.1 allows us to prove the following convergence 
theorem. 

Theorem 5.2 (Convergence as a 0). Let Uq , So , ug^ , ;Bg= ^ H^r\V. Let {\i,B) 
he the solution to (1.1) with v = 0, fi > 0, and with initial data (uo,/Bo). Let 
(u",S") be the solution to (1.2) with initial data (uo,Bq) . Choose an arbitrary 
T S (0, Tmax), where Tmax S (0, oo] is the maximal time for which a solution (u, B) 
exists and is unique. Suppose that Uq xit^ inV and Bq — > Bq in H . Then u" — > u 
m L°°([0,r],T/) andB"' ^B in L'^{[0,T],H) n L^{[0,T],V), asa^O. 

Proof. Under the hypotheses on the initial conditions. Theorem 5.1 holds, and 
thus there exists a time T > and a unique {u,B) e C{[0,T], H^{T'^) Ci V) x 
Ci[0,T],H^{T'^) n V) solving (1.1) (with ly = 0, fi > 0, (in particular, it holds 
that Tmax > 0). Thanks to Theorem 3.3, we know that there also exists a unique 
solution to the problem (1.2), namely (u",S") e Ci[0,T],V) x {L^{{0,T),V{A)) n 
C{[0,T],V)). Subtracting the equations corresponding to {u,B) and (u",S"), and 
recalling the fact that ^ = — A due to the periodic boundary conditions, we find 

(5.1a) -a2^Au" + 4(u"-u) 

dt dt 

= B[B" -B,B) + B{B'\ B" -B)- S(u" - u, u) - B(u", u" - u), 

(5.1b) -mA(6" -S) 

dt 

= S(S" - B, u) + B(S°, u° - u) - B(u" - u, S) - B(u", i3° - B), 

We now take the inner product of (5.1a) with u" — u and of (5.1b) with — B, 
and add the results. Using (2.3) and rearranging the terms, we find 

(5.2) (a^llu - u"f + |u - u"|2 + |i3 - S"p) + - Bf 

= -(S(u" - u, u), u" - u) - (S(S" - B, B),u" - u) 

- (B(u" -u,B),B" - B) - (B(6" - B, u), B" - B) + (Aut, u" - u) 

< |u-u"n|Vu||L.=c +2|S-B"|||VB|Uo.|u" -u| + |6-S"n|Vu|U=c 
+ (Aut,u" - u) 

< Ki{\u - u"|2 + |S - + (Aut, u" - u) , 

where Ki := Csupjg max{||Vu||icc, ||VS||l=°} < oo, and where we have used 
Young's inequality. It remains to estimate the last term on the right-hand side. 
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Using the fact that u satisfies (1.1) (with = 0), we have 
(5.3) {Aut,u"' 

= ( A[6 • VB - u • Vu] , u" - u) 

{AB ■ VB + 2{VB ■ V)VS + B ■ VAB, u" - u) 
- (Au • Vu + 2(Vu • V)Vu + u • VAu, u" - u) 
< Ca'\\B\\'/^BfJ^\\u - u"|| + Ca'mUl^ u"|| 

+ Ca^llull V^ilull-V^ilu - u"|| + Ca2||u||^.||u - u"|| 
</i2a^ + a^||u-u"f 
where X2:=Csup[o,T] max {||fi||V2||s||^/2,«2||g||2^^^|,„||i/2|,„||3/2^ 

Combining (5.2) with (5.3) yields 

^ (a^llu - u"||2 + |u - u"|2 + \B- S"|2) + 2n\\B" - B\\^ 
< K3{a^\\u - u"f + \u-u"f + \B- B"|2) + a^K2, 
where K3 ~ C max {1, 7^2}. Using Gronwall's inequality, we have 
a^\\u(t) - u"(i)||2 + |u(<) - u"(<)|2 + \B{t) - S"(OP 



< 00. 



2^ r||S"(s)-S(s)fe^^(*-^) ds 
Jo 



< C (a^lluo - <f + |uo - <P + \Bo BS\') + a'^il e^^*). 

Thus, if Uq uq m V and Bq ^ Z?o in iJ as a (in particular, if Uq = uq 
and B^ = Bq for all a > 0), then u" ^ u in L°°{[0,T],V) and ^ B in 

i°°([o,r],iJ) ni2([o,r],F) as a ^ 0. □ 

Theorem 5.3 (Blow-up criterion). With the same notation and assumptions of 
Theorem 5.2, and taking Ug = Uq and Bq ~ Bq for all a > 0, suppose that for 
some < 00, we have 

(5.4) sup limsupa2||u"(t)||2 > 0. 

te[0,T,) a^O 

Then the solutions to (1.1) with ly ~ 0, /i > become singular in the time interval 

[o,r,). 

Proof. Suppose for a contradiction that u and B are bounded in L°°([0, T], iJ'^(T^)n 
V), but that (5.4) holds. Taking the inner product of (1.2a) with u" and the inner 
product of (1.2b) with B" and adding the results, we find after using (2.3), 

«2||u"(t)||2 + |u"(i)|2 + |B"(i)p +^l f ||S"(s)f ds = a^u^r + |uoP + \Bo\'. 



Taking the lim sup as a — > 0^ gives 

(5.5) limsupa2||u(t)f + |u(i)p + |B(0|2 + /i f \\B{s)\\^ ds = \uo\^ + \Bo\ 

a^0+ Jo 
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However, thanks to Theorem 5.1, we have u,B E C{[0,T], CiV). Using Lions' 
Lemma (see, e.g., [58, p. 176] or [47, Corohary 7.3]) and (2.3) it is straight-forward 
to prove from (1.1) (with u = 0, fi > 0), that on [0,T], 

\u{t)\' + |B(t)p + fi f \\B{s)rds = |uo|2 + l^oP. 
Jo 

so that (5.5) contradicts (5.4). □ 

Remark 5.4. We note that that in the case fi,iy > 0, theorems similar to Theorems 
5.1, 5.2, and 5.3 can be proven with weaker assumptions on the initial data. In the 
case of Theorem 5.1, one can use, for example, Galerkin methods, and ideas similar 
to those of in the theory of the Navier-Stokes equations (see, e.g., [17, 58]). The 
analogous of Theorems 5.2 and 5.3 can be proven using nearly identical methods 
to those employed above. 
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